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Multiply periodic states appear in a wide variety of physical contexts, such as the Rayleigh-Benard 
convection, Faraday waves, liquid crystals, domain patterns in ferromagnetic films and skyrmion 
crystals recently observed in chiral magnets. Here we study a simple model of an anisotropic 
frustrated magnet and show that its zero-temperature phase diagram contains numerous multi-g 
states including the skyrmion crystal. We clarify the mechanism for stabilization of these states, 
discuss their multiferroic properties and formulate rules for finding new skyrmion materials. In 
addition to skyrmion crystal, we find stable isolated skyrmions with topological charge 1 and 2. 

Physics of isolated skyrmions in frustrated magnets is very rich. Their statical and dynamical 
properties are strongly affected by the new zero mode - skyrmion helicity. 

PACS numbers: 75.70.-i, 75.10.-b 75.30.Kz 


Introduction 

Magnetic skyrmion is a topological defect with a com¬ 
plex non-coplanar spin structure [T]. Predicted more 
than 20 years ago [2], skyrmions have been recently ob¬ 
served in conducting and insulating helimagnets under an 
applied magnetic field mil]. Magnetic skyrmion provides 
a physical realization of the idea that quantization of 
physical observables, such as electric and baryon charge, 
is rooted in topology 016 ]. The skyrmion spin texture in¬ 
duces one flux quantum of an effective magnetic field act¬ 
ing on spin-polarized electrons and magnons, which gives 
rise to topological Hall effects in charge and heat trans¬ 
port and at the same sets skyrmions into motion [THU]. 
Low critical currents needed to manipulate skyrmions 
opened a new active field of research - skyrmionics, which 
has a goal of developing skyrmion-based magnetic mem¬ 
ory and data processing devices pRUTB] . 

So far, skyrmions have only been found in a handful 
of materials, all with the non-centrosymmetric cubic B20 
structure, in which the relativistic Dzyaloshinskii-Moriya 
interaction turns a collinear spin state into a helical spi¬ 
ral. Under an applied magnetic field the spiral trans¬ 
forms into the skyrmion crystal (SkX) with a triangular 
magnetic superlattice. Further increase of the magnetic 
field results in a transition to the saturated ferromag¬ 
netic state, in which isolated skyrmions exist as stable 
topological defects [HI HE]. 

Skyrmions are close relatives of magnetic bubbles 
(cylindrical domains) and the transformation of the heli¬ 
cal spiral into SkX is analogous to the field-induced tran¬ 
sition between the stripe domain state and the bubble ar¬ 
ray taking place in thin ferromagnetic films. In the frame¬ 
work of Landau theory the bubble array is described by 
three coexisting spin modulations with the wave vectors, 
qi^q 2 and qs, which add to zero. This 3g-state is sta¬ 
bilized by a non-linear interaction between these mod¬ 
ulations and the uniform magnetization induced by an 


applied magnetic field m- In bulk chiral magnets SkX 
competes with the conical spiral state and an additional 
order-from-disorder mechanism - stabilization by thermal 
spin fluctuations - was invoked to explain why skyrmions 
are only observed at elevated temperatures [3 [20]. 

Practical applications of skyrmions crucially depend 
on finding new classes of skyrmion materials and new 
microscopic mechanisms for their stabilization. Recently, 
Okubo et al studied numerically an isotropic Heisenberg 
spin model on a (centrosymmetric) triangular lattice with 
competing spin interactions [21] . It shows a spiral ground 
state, which above a critical magnetic field transforms 
into the SkX. In addition, a state with two coexisting 
spirals was found in applied magnetic fields. These multi- 
q states are induced by thermal fluctuations and appear 
at nonzero temperatures. 

Here we discuss a new mechanism which lends stabil¬ 
ity to multiply periodic states even at zero temperature. 
We show that an uniaxial magnetic anisotropy strongly 
affects spin ordering in the frustrated triangular magnet: 
Tuning magnetic field and anisotropy at zero temper¬ 
ature, we find eight different phases, five of which are 
multi-^ states. A large part of the phase diagram is oc¬ 
cupied by the SkX. 

We clarify the nature of the phases found in numerical 
simulations and present results of analytical studies of in¬ 
stability of spiral states towards additional periodic mod¬ 
ulations. We obtain the spectrum of low-lying excitations 
in all ground states, which can be measured in magnetic 
resonance experiments, discuss their multiferroic proper¬ 
ties and the possibility to observe large magneto-electric 
responses in multi-g' states. 

We also explore very rich physics of isolated skyrmions 
in frustrated magnets. We find that their monopole 
and toroidal moment densities oscillate with the distance 
from the skyrmion centre, which gives rise to alternation 
of attraction and repulsion between skyrmions, cluster¬ 
ing of skyrmions in high magnetic fields and stability of 
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skyrmions with the topological charge 2. Dynamics of 
skyrmions in frustrated magnets is different from that in 
chiral magnets because of the additional collective degree 
of freedom - skyrmion helicity. The coupling between the 
helicity and the skyrmion center-of-mass motion leads to 
interesting dynamical effects. 

The model 

We consider classical spins, S^, of unit length on a 
triangular lattice in the x^-plane with ferromagnetic 
nearest-neighbour (NN) and antiferromagnetic next- 
nearest-neighbour (NNN) exchange interactions: 

E = -JiY,^i- Si + J2Y.^i' 

{i,i> ((i,i» 

( 1 ) 

i i 

where (i,j) and {{i^j)) denote pairs of NN and NNN 
spins, respectively, and Ji, J2 > 0 . The third and the 
fourth terms describe, respectively, the interaction with 
the magnetic field parallel to the z axis and the single-ion 
magnetic anisotropy. 

The critical ratio J2 / Ji for appearance of magnetic spi¬ 
rals and other modulated spin states is 1 / 3 . This can be 
understood by considering a modulation with the wave 
vector along the y axis parallel to a NNN bond (see inset 
in Fig. [^), in which case the model effectively reduces 
to a spin chain running in the y direction with the ferro¬ 
magnetic NN interaction, J[ = 2 (Ji — J2), and the anti¬ 
ferromagnetic NNN interaction, 7 ^ = 72. A spiral state 
in a chain appears for which corresponds to 

J2/ Ji > In what follows, the value of Ji is set to 1. 

Phase diagram 

Figure shows the zero-temperature phase diagram of 
the model in the (i^, h)-plane for J2 = 0 . 5 , which counts 
eight phases. Their spin configurations can be described 
in terms of the Fourier series, 

1 ^ 

= Ao + - ^ -h c.c.) -h higher harmonics, 

^ a=l 

( 2 ) 

where qi = g(0,1), q 2 = q{-^,-\) and qz = 
— |) are the wave vectors of the three fundamen- 
tal modulations minimizing the exchange energy (inset in 
Fig. [^), such that + g^2 + ^3 = 0 . In all states, except 
the state VII, the uniform spin component Aq is parallel 
to the applied magnetic field: Aq = (0,0, Ag). The am¬ 
plitudes of higher harmonics, i.e. modulations with other 
wave vectors from the triangular lattice in the reciprocal 



FIG. 1: Zero temperature phase diagram of the frus¬ 
trated triangular antiferromagnet. The eight phases are: 
(I) fully polarized ferromagnetic state, (II) conical spiral, (III) 
2g-state, (IV) 2g'-state, (V) vertical spiral, (VI) vertical spi¬ 
ral plus two in-plane sinusoidal modulations, (VII) flop state 
and (VIII) skyrmion crystal. Dashed line is the upper critical 
field, above which isolated skyrmions are unstable. The NNN 
exchange constant, J 2 = 0.5, magnetic anisotropy, A, and 
magnetic field, /i, are measured in units of the NN exchange 
constant, Ji — 1. 


space spanned on qi and gf2, are relatively small. The 
eight ground states of the model are: 

(I) The fully polarized ferromagnetic state with Ag = 1 
and Aq, = 0. 

(II) The conical spiral state, which occupies the region 
of A < 0 (easy plane anisotropy) and h < h^{K) = 

This is a single-g' state with only one 
non-zero describing a circular spiral in the horizontal 
plane, e.g., A\ = ±iAf and Af = 0, ± sign corresponding 
to two opposite directions of spin rotation in the spiral. 

(III) The 2 g-state with |Ai| > IA2I > \As\ (see Fig.|^). 
For A > 0 (easy axis anisotropy) the conical state ac¬ 
quires the second spiral modulation in the horizontal 
plane, e.g., with the wave vector q2 and the opposite 
sense of spin rotation compared to that in the first spi¬ 
ral: A2 = TiA2,A2 = 0 . This ‘ 2 g-state’ has also a 
small sinusoidal spin modulation parallel to the 2: axis 
with the wave vector q^: A3 = (0,0, Ag). As explained 
in the next section, the sinusoidal modulation stabilizes 
the ‘ 2 g-state’ state. For small positive A, A2 oc \/A, i.e. 
A = 0 is a critical line. 

(IV) The 2 g'-state, which is the same as the state III but 
with |Ai| = IA2I (see Fig. [^). As A increases, the am¬ 
plitude of the second spiral grows until it becomes equal 
to that of the first spiral. This transition is marked by 
dotted line in Fig. 

(V) The spiral in a vertical plane, e.g., Ai = 
{i cos(x)A||, i sin(x)A||, A^) and Ay « A_l, where the an¬ 
gle X describes the rotation of the spiral plane around 
the 2: axis. This is the ground state for low applied fields 
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and K > ^. The vertical spiral is not perfectly circular 
and has higher harmonics. 

(VI) The spiral in a vertical plane with two sinusoidal 

modulations in the direction normal to the spiral plane, 
e.g., the spiral in the yz plane with the wave vector 
qi and Ai and the sinusoidal modula¬ 

tions with the wave vectors q 2 and qs along the x di¬ 
rection: A 2 = (^2,0,0) and A3 = (—A2,0,0). The sinu¬ 
soidal modulations appear above a critical magnetic field, 
hi{K)^ and the transition between the states V and VI 
is of second order. 

(VII) The ‘flop state’ (Fig. [^). As magnetic field in¬ 
creases further, the spiral plane changes from vertical to 
horizontal. This flop transition is not instantaneous, but 
occurs in a narrow field interval, h 2 {K) < h < h^{K)^ in 
which the spiral plane as well as the direction of the sinu¬ 
soidal modulations rotate. The orientation of the spiral 
plane changes abruptly at ^2 and at h^{K) the system 
undergoes a first order transition to the 2g-state. The 
flop state is the only state, in which the uniform magne¬ 
tization Aq has a non-zero component parallel to the xy 
plane. 

(VIII) The triangular skyrmion (or antiskyrmion) crys¬ 
tal state (Fig. [^) with equal amplitudes of the three 
fundamental modulations: 

Aa = ^e*'^“(-zsin(xa),icos(x„),l). (3) 

These modulations are spirals in three vertical planes 
rotated with respect to each other by ±120°: = 

X ± v{a — 1)^ {a = 1,2,3), where v = ±1 is the 
skyrmion vorticity describing the sense of rotation of 
the in-plane components of spins along a contour encir¬ 
cling the skyrmion centre [1]. The sum of the phases 
01 ± 02 ± 03 is either 0 or tt and the ^-component of spin 
in the centre of each skyrmion is = cos(0i ± 02 ± 03)- 
The topological charge, Q, of skyrmions forming the crys¬ 
tal is given hy Q = v cos(0i ±02 ±03). The spin direction 
in the skyrmion centre is opposite to that of the applied 
magnetic field. Because of the arbitrary sign of vorticity, 
the topological charge of skyrmions can be both ±1 and 
— 1 for any direction of the field. This degree of freedom 
does not exist for skyrmions in chiral magnets and mag¬ 
netic bubbles [21]. The SkX occupies a significant part 
of the phase diagram and all its boundaries are lines of 
first-order transitions. 


Instabilities of spiral states 

To clarify the origin of multi-g states, we discuss a 
continuum model of the frustrated triangular magnet, 
applicable when the modulation period is much larger 
than the lattice constant (g <C 1): 

E = Jd^xi^^SDS+^iSn^-hS^y (4) 




X 



FIG. 2: Multi-g states. In-plane components (arrows) 

and out-of-plane components (color) of spins in four selected 
states from the phase diagram Fig. (a) 2g-state (state 
III), (b) 2g'-state (state IV), (c) Flop state (state IV) and 
(d) Skyrmion crystal (state VIII). The inset shows the wave 
vectors of the three fundamental modulations in the multi-g 
states. 


where the first term is the exchange energy, the second 
term is magnetic anisotropy and the last term is the Zee- 
man energy. The operator D in the momentum repre¬ 
sentation can be expanded in powers of the wave vector 

Q‘ 


Dq — do ± d2q^ ± d4g^ ± dgg^ 

+ d'e {qt - 'i-^qUy + - 9 ®) + • • • • ( 5 ) 

The anisotropy in the reciprocal space first appears in 
the sixth order of the expansion, resulting in 6 minima 
of Dq 8 it q = ±gri, ±gr2 5 i^3- It is convenient to add a 
constant to energy to make Dq^ = 0. 

The six-fold degeneracy of the exchange energy min¬ 
imum is the source of complexity. Another important 
factor is the relative rigidity of single-spiral states, which 
for K^h ^ 0 makes the states with several coexist¬ 
ing modulations, spiral or sinusoidal, energetically more 
favourable, even though the multi-g states have a larger 
exchange energy due to higher harmonics induced by the 
local constraint S‘^{x) = 1. 

To prove this point, we consider the case of weak 
anisotropy, \K\ 1, and magnetic fields just below 

h^{K) = Dq — K, at which the transition to the satu¬ 
rated ferromagnetic state occurs. Then the in-plane com¬ 
ponent, m, of spin S = m A is small and the 

energy can be expanded in powers of m: 

epM + ^’’^0 + dh)m 

+ lm^{D + 6h- K)m^ , 

8 



(6) 
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where 5 fm is the energy density of the saturated ferro¬ 
magnetic state and 6h = h — h^{K) < 0 . For K < 0^ 
the ground state of the model is the conical spiral (state 
II): mo = A {x cos{qi'X a) ± y sm{qi'X a)) with 


A = '\J arbitrary phase a. For K > 0, 
the conical spiral state becomes unstable against forma¬ 
tion of an additional spin modulation. A small deviation 
from the conical spiral, k = m — mo, changes energy by 
SE ^ I f£x jkBk + {k-mo){D - K){k-mo)\ . An ad¬ 
ditional modulation satisfying 


I Dk = 0 , 

\ D{k’mo) = 0 , 


( 7 ) 


changes energy hy SE = Jd‘^x{kmo)^ < 0 , for A" > 0 . 
The solution of Eqs. 0 is the second x^-spiral with the 
wave vector q2 (or q^) and the opposite direction of spin 
rotation: k = B {x cos{q2' x/S) ^ y sm{q2 'x/S)). 
The two-spiral state III has a lower energy because of 
the small sinusoidal modulation of the 2:-component of 
spin with the wave vector q^^ SS^ = —{k- mo) = 
—AB cos{qs ' X — a — necessary to preserve the spin 
length. This sinusoidal modulation also has the minimal 
exchange energy, but its anisotropy energy for iF > 0 is 
lower than that of the xy spirals, which makes the conical 
state unstable. 

Next we discuss the instability of the vertical spiral 
(state V), 


S = ±y sin 0 -h i cos 0 , 


(8) 


for h,K 1 . The magnetic field and easy axis parallel 
to the spiral plane deform the spiral: (j) ~ qi'X — asm{qi • 
x) — 6 sin { 2 qi • x) with a « D0+D2 — ^ ^ 

responding to 

^ ^^cos{qi’x)-^cos{ 2 qi’x)-^cos{ 3 qi’x). ( 9 ) 

Importantly, the uniform magnetization equals the am¬ 
plitude of the energetically costly second harmonic. This 
makes the magnetic susceptibility of the ^2:-spiral low 
and leads to appearance of an x-component of spin 
above a critical magnetic field. Expanding the spin 
energy in powers of around Eq. (|^, we obtain 
sm{(t))DS^ — S^D sin(0) = 0, the solution of which is the 
sum of two sinusoidal spin modulations with the wave 
vectors q2 and q^ plus relatively small higher harmonics: 

^ C [cos(gr2 • X + 7) — cos(gr3 • cc — 7)] , 7 being an ar¬ 
bitrary phase. The second-order transition to the non- 
coplanar 3 g-state VI occurs at hi{K) ^ , 

when the lowering of the Zeeman energy due to the larger 
magnetization of the 3 g-state compensates the increase 
of the anisotropy energy. 

Similarly, the SkX owns its stability partly to the low 
anisotropy energy (as it is composed of three vertical spi¬ 
rals) and partly to its large magnetic moment. The latter 



FIG. 3: Helicity reversals, (a) False colour plot of the 
toroidal moment density in the skyrmion with the helicity 
angle y = 7r/2. (b) Fan-like oscillations of the out-of-plane 
spin component . The skyrmion centre is located at x = 0. 


is the consequence of the constraint = 1, which is im¬ 
possible to maintain without adding to the 3 g-state a 
large uniform magnetization. 

The discussion of instabilities of spiral states in the 
continuum model of the frustrated triangular magnet 
also shows that the phase diagram Fig. is generic, a 
minor difference being that in the discrete model with 
J2I J\ — 1/2 the wave vectors of the three fundamen¬ 
tal modulations are commensurate with the lattice and 
constant throughout the phase diagram, while in the con¬ 
tinuum model g is a function of K and h. 


Helicity reversals and interaction between skyrmions 

Next we discuss properties of isolated skyrmions. The 
asymptotic of the in-plane spin components at a large 
distance r from the skyrmion centre is given by 

( u(r) ( ^ ^ flO) 

y J 7 V sin <p cos (p J y sin X J ’ 

where (p is the azimuthal angle in the xy plane and v and 
X are, respectively, the skyrmion vorticity and helicity 
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angle introduced in Eq. © (for simplicity, we consider 
large skyrmions and replace the discrete rotational sym¬ 
metry of the triangular lattice by the continuous one). 

In absence of in-plane magnetic anisotropy, the helicity 
angle y is arbitrary, y = 0 , tt corresponding to skyrmions 
carrying a monopole moment, A oc X]j 
y = ±7r/2 corresponding to skyrmions with a toroidal 
moment, oc x SjY | 22 ]- Arbitrary y is a zero 

mode, which skyrmions in chiral magnets and magnetic 
bubbles do not have. 

In addition, the toroidal and monopole moment den¬ 
sities show periodic sign reversals as the distance r from 
the center increases (see Fig.[^, similar to the helicity re¬ 
versals observed in magnetic bubbles in a ferrimagnetic 
hexaferrite, where y can be both +7r/2 and —7r/2 [23] . 
We find it more convenient to keep the helicity angle fixed 
and associate these oscillations with the sign changes of 


u{r) in Eq. (10), 


/N A 
uir) ^ V 

Vr 


cos 



( 11 ) 


where g ~ y ^ is the length of the minimal-energy wave 

vector and The sign reversals of u{r) corre¬ 

spond to fan-like oscillations of spins around the z axis, 
whose amplitude decreases exponentially with r for k > q 
or, equivalently, h > h^{K). This condition ensures that 
the skyrmion energy is finite and, hence, the lower criti¬ 
cal field for stability of isolated skyrmions coincides with 
the line separating the fully polarized state from the con¬ 
ical spiral state, for K <0, and the 2^'-state, for K >0. 
There is also an upper critical field (dashed line in Fig.[^, 
above which skyrmions collapse. 

The fan oscillations of spins in skyrmion give rise to 
sign changes of the skyrmion-skyrmion interaction U (^ 12 ) 
(Fig. [^). This interaction depends on the topological 
charges of skyrmions. For brevity, we call a skyrmion 
with positive vorticity and Q = — 1 a skyrmion, while a 
skyrmion with v = —1 and Q = -hi is called an anti- 
skyrmion. Two skyrmions or two antiskyrmions attract 
each other at distances of the order of the skyrmion di¬ 
ameter, the maximal reduction of energy being ^ 10% 
of the skyrmion energy. Because of the attraction, the 
skyrmion crystal phase extends to h > h^{K)^ as the 
energy per skyrmion in the crystal can be negative even 
when the energy of an isolated skyrmion is positive. Iso¬ 
lated skyrmions tend to form clusters with a short-range 
crystal order, resembling the clustering of vortices in 
1.5-type superconductors [24]. The oscillating interac¬ 
tions between skyrmions in the frustrated magnet are in 
sharp contrast with skyrmion-skyrmion interactions in 
non-centrosymmetric magnets, which do not show fan¬ 
like oscillations and repel each other at all distances [T8| . 

Interactions between skyrmions also depend on their 
helicities: t/i 2 (^) ~ V{r) + W{r)cos{xi — X 2 ) (see 


Fig. 1^), which couples the helicity dynamics to the 
translational motion of skyrmions. Another difference 
from chiral magnets is the existence of the locally sta¬ 
ble skyrmion with topological charge Q = ±2, shown 
in Fig. [^, which can be considered as a tightly bound 
state of two skyrmions with opposite helicities, similar 
to the magnetic bubble with Q = 2 observed in thin 
films of a CMR manganite [25] . The energy of the Q = 2 
skyrmion is, however, larger than the energy of two Q = 1 
skyrmions. 

The helicity dependence of the skyrmion-antiskyrmion 
interaction is more complex: it is approximately propor¬ 
tional to cos(2(asa + Xs — Xa), whcrc asa is the angle 
between the line connecting skyrmion with antiskyrmion 
and the x axis, with respect to which the skyrmion and 
anstiskyrmion helicity angles, Xs Xa, are defined. 
Figure]^ shows skyrmion-antiskyrmion interactions for 
Xs = Xa and aga = 0. The minimum of the skyrmion- 
antiskyrmion potential occurs at larger distances and is 
more shallow than that of the skyrmion-skyrmion po¬ 
tential. One can form a metastable rectangular crystal 
with alternating columns of skyrmions and antiskyrmions 
(Fig.|^). Such a crystal with zero topological charge has 
a higher energy than a purely skyrmion (or purely anti¬ 
skyrmion) crystal, which is important for observation of 
Topological Hall Effects in frustrated magnets. 


Magnetoelectric coupling 

All modulated states in the phase diagram Fig.l^break 
inversion symmetry and can induce an electric polariza¬ 
tion, P, in a magnetic insulator. The two possible direc¬ 
tions of spin rotation in spirals correspond to two oppo¬ 
site directions of P [26] . 

Consider, for example, a quasi-two-dimensional mag¬ 
net with the ‘eclipsed’ triangular spin layers stacked on 
top of each other m and the symmetries of the crys¬ 
tal lattice being inversion /, three-fold axis ‘^z and the 
two-fold axis 2^ (e.g. space groups P3mi, R3m). This 
allows for three independent Lifshitz invariants describ¬ 
ing the magnetically-induced electric polarization, which 
in momentum representation have the form: 

Hue = [giE^Sl^iq^+ qySy) 

Q 

+g2{q^Ey - qyE^sigSy ( 12 ) 

+938!^ {{q^E^ - qyEy)sy + {q^Ey + qyE^)Sl) ]. 

The polarization, P = — , described by the first two 

terms lies in the spiral plane and is orthogonal to the spi¬ 
ral wave vector, as is the case for P induced by a cycloidal 
spiral through the inverse Dzyaloshinskii-Moriya mecha¬ 
nism [28ll3Q] . The third term describes other symmetry- 
allowed directions of polarization, in particular, P par¬ 
allel to the wave vector of a helical spiral m- 
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(a) 



(b) (c) 



FIG. 4: Unusual features of skyrmions in frustrated 
magnets, (a) The interaction energy, U 12 , vs the distance ri 2 
between two skyrmions (black lines with rectangular mark¬ 
ers), and skyrmion and antiskyrmion (red lines with circular 
markers). Solid (dashed) lines show the interaction for equal 
(opposite) helicities of the topological defects. Skyrmion posi¬ 
tions and helicities were fixed by the condition = — 1 at the 
skyrmion centre and Oi = 0 ai the six nearest-neighbour sites, 
(b) Metastable skyrmion with topological charge Q — 2. (c) 
Metastable skyrmion-antiskyrmion crystal with a rectangular 
lattice. In panels (b) and (c) color indicates the topological 
charge density, pg, while arrows show the in-plane compo¬ 
nents of spins. 


Figure shows magnetic field dependence of the elec¬ 
tric polarization in arbitrary units, for K = 0.04, p 2 = 9\ 
and p 3 = 0. The states with spirals in a vertical plane, 
including the SkX, induce polarization in the ^-direction 
(black line), while the states with spirals in the horizontal 
plane (conical spiral and the two-spiral states) induce an 
in-plane polarization (red line). In the flop phase the spi¬ 
ral plane rotates in an applied magnetic field and so does 
P. The spontaneous in-plane magnetization, present in 
this phase (see Fig. Hf) allows for effective control of the 
polarization direction by magnetic field and vice versa. 



FIG. 5: Multiferroic and dynamical properties of 

multi-g states, (a) Magnetic field dependence at iF = 0.04 
of (a) the average spin vector {S), (b) the electric polarization 
vector P and (c) the imaginary part of the in-plane magnetic 
susceptibility, x"{^) (false color plot). Black and red lines in 
panels (a) and (b) are the out-of-plane and in-plane compo¬ 
nents of th e ve ctors. The electric polarization was calculated 
using Eq. |T^ for gi = ^2 and ^3 = 0. The inset in panel (c) 
shows magnetic modes in the 2g-state near the transition to 
the flop state (see text). 


The orientations of spiral planes and skyrmion helicity 
in realistic materials will be affected by the sixth-order 
in-plane magnetic anisotropy, which we neglected. The 
interplay between this weak anisotropy and the inter¬ 
action of the electric polarization induced by the spiral 
with an applied electric field can lead to giant magneto¬ 
electric effects, such as the electrically-induced rotation 
of the skyrmion helicity. 


Magnetic resonance 

We also studied collective spin dynamics for the states 
in the phase diagram Fig. [^induced by the ac magnetic 
field, sin(u;t), parallel to the xy plane. Figure]^ shows 
the false color plot of the imaginary part of the in-plane 
dynamical magnetic susceptibility calculated along the 
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K = 0.04 line in the phase diagram (as in panels (a) 
and (b)) by solving Landau-Lifshitz-Gilbert equation (see 
Methods). The phase transitions between different states 
are reffected in sudden changes in the magnetic resonance 
spectrum. 

While the SkX in chiral magnets shows two resonances: 
one with the counterclockwise and another with the anti- 
counterclockwise rotation of the skyrmion centres [38] , in 
the frustrated magnet we find only one mode, in which 
skyrmions rotate counterclockwise. Figure shows cir¬ 
cular trajectories traced by the skyrmion centres for vari¬ 
ous amplitudes of the ac field and the resonant frequency 
cj = 0.224Jifi-h 

The field-induced spin dynamics in frustrated magnets 
is, however, more complex than in chiral magnets, be¬ 
cause of the additional zero mode - helicity. This can be 
seen from the fact that the time dependence of an in¬ 
plane component of spin at a chosen site (see Fig. I#) is 
a superposition of periodic oscillations with two different 
periods. The shorter period corresponds to the rotation 
of the skyrmion centre, while the longer one results from 
a monotonous increase of the helicity angle, as evidenced 
by the series of snapshots shown in Fig. (see also the 
movie in Supplementary Material). The helicity dynam¬ 
ics is induced by the rotation of skyrmion centres through 
a non-linear coupling between the two modes, as can be 
seen from the non-linear dependence of the frequency of 
helicity rotations, uJheh on tko amplitude of the ac mag¬ 
netic field (Fig. The coupling between the helicity 

and translational dynamics must have important conse¬ 
quences for the current-induced motion of skyrmions in 
confined nanostructures. 

Another remarkable feature of the magnetic resonance 
spectrum is the presence of the zero-frequency peak in 
the fiop phase (Fig. [^) related to the non-zero spon¬ 
taneous in-plane magnetization, which in absence of the 
in-plane magnetic anisotropy has an arbitrary direction. 
Surprisingly, a low-energy mode is also present in the 2q- 
state near the first-order transition to the fiop phase, as 
shown in the inset to Fig. (this transition does not 
actually occur at K = 0.04 because of the intervening 
SkX phase, which was artificially removed to obtain the 
susceptibility shown in the inset). 


Conclusions 

In conclusion, we studied a new class of skyrmion mate¬ 
rials - frustrated magnets - and clarified the mechanism 
for stabilization of various multi-^ states which does not 
involve thermal spin fiuctnations. The crucial role in this 
mechanism is played by magnetic anisotropy, which to¬ 
gether with an applied magnetic field leads to instability 
of spiral states. An interesting question is why the 2q 
and skyrmion crystal states are favoured by both the en- 
tropic mechanism in absence of anisotropy m and the 



FIG. 6: Coupled dynamics of skyrmion helicity and 
centre-of-mass. (a) Circular trajectories spanned by the 
centres of rotating skyrmions in the in-plane ac magnetic field 
for various amplitudes of the ac magnetic field huj. (b) Fre¬ 
quency of the helicity rotation, UheU rneasured in units of the 
frequency of the ac field uj = 0.224Ji/i“^ vs huj. (c) Time evo¬ 
lution of the x-component of spin at a chosen site, marked by 
a circle in the series of snapshots (panel d), which also shows 
current positions of the skyrmion centre (red point) and its 
trajectory for hco = 0.01. This calculation was performed for 
h = 0.2 and K = 0.04. 


magnetic anisotropy at zero temperature? The common 
origin seems to be the soft magnetic excitation spectrum 
of multi-g states, which on the one hand lowers free en¬ 
ergy of these states by increasing their entropy and, on 
the other hand, lowers their energy by enhancing suscep¬ 
tibility to magnetic field and anisotropy. Soft magnetic 
modes of multi-g states, which can be excited by both 
electric and magnetic fields, lead to large magnetoelec¬ 
tric responses. 

Finally, we list a few requirements for stabilization 
of skyrmions in frustrated magnets following from our 
study, which can be used as guidelines in the search for 
new materials hosting skyrmions: (1) spiral spin order¬ 
ing in zero magnetic field resulting from competing ex¬ 
change interactions; (2) three-fold or six-fold symmetry 
axis, which ensures that the spirals with the three wave 
vectors, qi, q 2 and qs forming the skyrmion crystal, have 
the same energy; although cubic lattice symmetry does 
not favour the three vectors lying in the same plane, 
it results in a weak dependence of the exchange energy 





on gr, which may be sufficient; (3) easy axis magnetic 
anisotropy axis stabilizing multi-g states. Although these 
requirements are not very demanding and some of them 
are met in NiGa 2 S 4 [32] [33], (a-NaFe02 [34], dihalides 
NiBr 2 and Fea^Nii_a,Br 2 [354137] . it is a challenge to find 
a material which satisfies them all. 


Methods 

The energy 0 has been minimized using the iterative 
simulated annealing procedure and a single-step Monte- 
Carlo dynamics with the Metropolis algorithm (for de¬ 
tails see Ref. [H]). We imposed periodic boundary con¬ 
ditions and performed simulations for lattices of different 
sizes to check the stability of the numerical routine. 

The dynamical magnetic susceptibility for all multi-g 
states and the time evolution of spins in the oscillating 
magnetic field have been obtained by solving the Landau- 
Lifshitz-Gilbert (LLG) equation using the fourth-order 
Runge-Kutta method. The dimensionless damping pa¬ 
rameter a = 0.01. To obtain initial spin configurations, 
we used the annealing technique supplemented by a fur¬ 
ther LLG relaxation for = 0. When the convergence 
of the spin configurations was reached, we switched on 
a periodic magnetic field in the plane xy to study the 
spin dynamics (Fig. or applied a magnetic field pulse 
at t = 0 and monitored the response, to obtain the mag¬ 
netic susceptibility (Fig. [^). 

The skyrmion-skyrmion and skyrmion-antiskyrmion 
interaction potentials shown in Fig.[^ were calculated by 
minimizing the spin energy with the constraint, = —1, 
imposed at the centres of the topological defects. To 
calculate the helicity dependence of the interactions, we 
constrained the helicity angle at six sites neighbouring to 
the skyrmion centre and imposed = —0.427 at these 
sites, which holds for an isolated skyrmion. 
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